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\Q Abstract 






The two-dimensional one-component plasma at the special coupling /3 = 2 is known to 
be exactly solvable, for its free energy and all of its correlations, on a variety of surfaces 
and with various boundary conditions. Here we study this system confined to a spherical 
annulus with soft wall boundary conditions, paying special attention to the resulting 
asymptotic forms from the viewpoint of expected general properties of the two-dimensional 
plasma. Our study is motivated by the realization of the Boltzmann factor for the plasma 
system with /3 = 2, after stereographic projection from the sphere to the complex plane, 

^^ by a certain random matrix ensemble constructed out of complex Gaussian and Haar 

(^ distributed unitary matrices. 

(N 

^ 1 Introduction 

o 

,-H The two-dimensional one-component plasma is an equilibrium statistical mechanical system 

^ consisting of A^ mobile particles, each of charge +1, and a smeared out neutralizing background. 

>> The particles are confined to a two-dimensional surface, and the charge densities (both point 

S^ and continuous) interact through the solution of the two-dimensional Poisson equation on the 

^ surface. 

Although it is defined as a classical system, the two-dimensional one-component plasma in 
the case that the surface is of constant curvature is also known in quantum many body physics. 
This is due to its relevance to the fractional quantum Hall effect. Thus it turns out that the 
Boltzmann factor for the plasma system at inverse temperature (3 = 2v, v an odd integer, 
is equal to the absolute value squared of the Laughlin trial wave function for the fractional 
quantum Hall effect at filling fraction \jv [351 [2H1 E]. In the case v = \ and thus /3 = 2 the 
corresponding trial wave function is in fact the exact wave function for non-interacting spinless 
fermions with constant perpendicular magnetic field. 

It has been known for some time that there is also an analogy between the two-dimensional 
one-component plasma confined to a disk in the plane, and the complex Ginibre random matrix 
ensemble [1]. The latter is specified as the eigenvalue probability density function (PDF) 



for N X N complex Gaussian matrices, where each element is independently distributed as a 
standard complex Gaussian. In terms of the notation Zj = Xj + ii/j, Xj,yj G M, it has the 
explicit form 

N 
1=1 l<j<k<N 



up to proportionality. If the extra condition that \zi\ < yN is imposed, then (1.1) is propor- 
tional to the Boltzmann factor for the one-component plasma at coupling /3 = 2, confined to 
a disk of radius R = \/N. Without this constraint, the eigenvalues are to leading order still 
confined to a disk in this radius (an example of the circular law p3| El [211 EH])- 

More recently, an analogy between two other random matrix ensembles and the one-component 
plasma confined to the other homogeneous constant curvature two-dimensional surfaces — 
namely the sphere and pseudosphere — has been specified. Thus in [M] it was shown that the 
eigenvalue PDF for random matrices A~^B, where A and B are independent complex Gini- 
bre matrices, coincides with the Boltzmann factor for the one-component plasma at /3 = 2 
on the sphere, after a stereographic projection of the latter. And in [19j it was shown that 
the eigenvalue PDF of truncations of unitary random matrices |12] has the same form as the 
Boltzmann factor for the one-component plasma on the pseudosphere at /3 = 2, after projection 
of the latter onto the Poincare disk. These examples of the one-component plasma had earlier 
been identified as exactly solvable two-dimensional statistical mechanical systems [6], [161 [32] . 
As an aside we mention that the one-component plasma confined to a surface of non-constant 
curvature — Flamm's paraboloid which occurs as the spatial part of the Schwarzschild metric 
from general relativity in two- dimensions — has recently been shown to also be exactly solvable 
at /3 = 2 [H], although as yet no random matrix analogy has been found. 

A topic of much current interest in random matrix theory is ensembles formed from the 
product f/y^/^, where [/ is a unitary random matrix and Y is positive definite [2Zl[ini[2Sl[S]- The 
motivation behind our work is to relate, for a particular class of random matrices Y generalizing 
the ensemble A~^B, an eigenvalue PDF obtained in this setting to the two-dimensional one- 
component plasma at /3 = 2 confined to a spherical annulus. The system is exactly solvable, 
being an example of a determinantal point process. Moreover, we will see that the asymptotic 
forms of the partition function, one and two point correlations, and the distribution of a general 
axially symmetric linear statistic all illustrate physical properties of the point process which 
are expected to hold for the plasma system in the same geometry but with /3 > |13j. 

In Section 2 the Boltzmann factor for the one-component plasma confined to a spherical 
annulus is calculated, as is its form upon a stereographic projection. In the case (3 = 2, and 
with the area of the spherical caps outside the spherical annulus certain rational fractions of 
the area of the sphere, a realization of the projected functional form of the Boltzmann factor 
as the eigenvalue PDF of a random matrix ensemble is given in Section 3. In Sections 4 and 
5 the plasma system at /3 = 2 is studied as an exactly solvable statistical mechanical model, 
and the corresponding large N asymptotic forms are computed and used according to the final 
sentence of the above paragraph. 



2 The plasma system 

Consider a sphere 5* of radius R, and let < ^ < vr refer to the usual azimuthal angle, and 
< < 27r refer to the polar angle. For two points {6,4>) and {0',<f)') on the sphere, let a 
refer to their relative angle when considered as vectors in M^. We know that the solution of the 
charge neutral Poisson equation 



2i?2 



(the sphere being a compact surface, charge neutrality is a necessary condition for existence of 
a solution), where Ss{{0, 0), (6*', 0')) is the delta function on the sphere, is then given by [H] 

<l>((^,0),(^',0')) = -log(2i?sin(a/2)). (2.1) 

Introducing the Cayley-Klein parameters, 

u := cos{e/2)e''^/^, v := -2sin(e/2)e-^^/2 ^2.2) 



we know (see e.g. P^ eq. (15.108)]) that (2.1) can be rewritten 



$((^,0),(^',0')) = -\og{2R\u'v-uv'\). (2.3) 

Let us mark two circles on the sphere corresponding to the azimuthal angles 9q and n — 6q, 
with < 9q < IT — 6g < 71. The surface of the sphere between these circles defines a spherical 
annulus. Let ^[o,6»q] denote the area of the spherical cap above 6q and thus including the north 
pole, and let A[7r-6i,,7r] denote the area of the spherical cap below n — 9q and thus including the 
south pole. We parametrize 6q and 6g by introducing Q and q such that 

Ao^Oq] ^ Q A^-e,M ^ g (2 A) 

47ri?2 i + g + g' 47ri?2 1 + g + Q- ^ ■ ) 

The plasma is specified by requiring that within the spherical annulus there be N mobile 
particles of charge +1 and a uniform neutralizing background. Both the discrete and continuous 



charges are to interact via the potential (2.1 ). It follows from (2.4) that the area of the annulus 



^[eQ,7r-6»,] is such that the uniform neutralizing background charge density is equal to 



47ri?2 



:i + Q + q)=:-p,. (2.5) 



We would like to compute the potential energy V{9') of the interaction of a particle at 
{6', (f)') in the spherical annulus, and the neutralizing background. For this purpose we extend 
the background to have uniform charge density — p^ throughout the sphere. To compensate, we 
must impose a uniform charge density pb in the spherical caps above 9q and below n — 6q. We 
can now proceed to compute the sought potential. Throughout we will ignore the 2R factor 



in the logarithm of (2.3): by charge neutrality, we can check that it must contribute a factor 



(2i?)^/3/2 to the Boltzmann factor. 



Proposition 2.1 We have that 



nf nt 

V{e') = Cn-NQ log sin -^-Nq log cos -, 



where 



l I l + (y + g 2 



'^'l + Q).og ^ + « 



l + Q + g 



(2.6) 



(2.7) 



Proof. The potential of the interaction of a particle with the uniform background covering 
all the sphere is independent of the location of the particle. Choosing this location to be the 



north pole, we see from (2.3) and the fact that on the surface of a sphere dS = R sin9 



that the corresponding potential energy is 

PbR'^ / sin 61 (log sin -j(i6' / 
Using the integral evaluation 



(2.8) 



/ sin X i log sin — ) dx = i —1 + 2 log 



sm - sm 



(2.9) 



with t = TT we see that (2.8) simplifies to 



--{l + Q + q). 



(2.10) 



Consider next the potential between a particle and the charge density pb in the spherical 
cap above 6q. This is equal to 



'Q 

-PbR' I de sin e 

'o Jo 



2lT 



log \uv — uv'\. 



(2.11) 



Simple manipulation gives 



a ni 

log \u'v — uv'\ = log COS - + log sin — + log 



1 _ tang/2 ^_,(^„^,) 
tan0V2 



(2.12) 



Note that the ratio of tan functions has magnitude less than one. Substituting into (2.11) 



vanishes, and hence (2.11) reduces to 



this latter fact implies the third term in (2.12) does not contribute since the integral over 



PbR^ (27r) / sin 6^ f log cos - + log sin — ) d9. 
Jo ^ 2 2 / 



(2.13) 



This simplifies by noting from (2.4) that 



R^(2n) [ '^ sine de = UnR^) ^ , 

Jo 'l + Q + q' 



(2.14) 



while use of (2.9) shows that 



sm U loff cos - du = — sm cos — losr cos — 

^2 2 2^2 



Q 



l+g , 1+g 
log 



1 + Q + q 1 + Q + q "1 + Q + g' 



(2.15) 



where the second equality follows by making use of (2.14). Substituting (2.14) and (2.15) 



in (2.13) we conclude that the potential between a particle and the charge density pb in the 



spherical cap above 9q is equal to 



6' N 



NQ log sin - + -(^Q + {l + q)\og- 



l + q 



+ Q + q 



(2.16) 



Replacing qhj Q and 6' hj n — 6' gives that the potential between a particle and the charge 
density pb in the spherical cap below n — 9q is equal to 



f) N / 1 

-iVglogcos- + —(q + {l + Q)\og — 



+ Q + q 



Adding together (2.10), (2.16) and (2.17) gives (2.6). 



(2.17) 
D 



Note that an equivalent viewpoint on the result (2.6) is that the potential 

9 9 

—NQ log sin - — Nq log cos - 

results from charges NQ and Nq at the north and south poles respectively. With aj denoting 
the angle between a point (^, 0) on the sphere, and another point (6'j, 0j), a related question is 
to seek the background charge density which gives rise to the potential 



i=i 

In a disk geometry, the analogous question has recently been addressed in [1]. 

We turn our attention next to the computation of the potential for the interaction of the 
background with itself. 

Proposition 2.2 The background-background potential is equal to 



N^ N^ , ,, 1 + g 

t2 



1 + Q) log 



l + Q + q 



N^ / 1 

+ —[-{Q + q) + Q\og—-—^+Q{l + Q)\og{l + Q)\og{l + Q)-QHogq 



glog 



l + Q + q 



Q + q 

+ g(l + q) log(l + q) log(l + q) - q^ log Q 



(2.18) 



Proof. The background-background potential is given in terms of the particle background 
potential V{9) according to 



-Pbi2nR^ 



sm9V{e)de. 



(2.19) 



Substituting (2.6) and performing the first of the resulting integrals gives 



N N 

--C^ + j{l + Q + q) 



e e^ 

sin 9 ( NQ log sin - + Nq log cos - ) dO. 



The integrals can be performed using (2.9) and further reduced as in the second equality of 
( |2.15D , with the result being ( |2.18D . D 

The total potential energy U of the plasma system consists of the particle-particle, particle- 
background, and background-background interactions. It therefore follows from (2.3), (2.6), 
(2.18) and the remark above Proposition 2.1 that the Boltzmann factor e~^^ for the plasma 

(2.20) 



system is equal to 



1 x^'/^/a 



{-k) 



N 



-/3X^-Q|^^|/3Q^|„^|/3<?iV -Q 



\UkVj -UjVkf, 



1=1 



^<j<k<N 



where 



N'^ / 1 

Kn : = — - (1 + Q + g) + 2(1 + Q + g) log + (1 + qf log(l + q) 

4 V l + U + a 



+ Q + q 
+ {l + Qf log(l + Q)-Q' log q - q' log q) . 



(2.21) 

By construction the particles are restricted to the spherical annulus. However, as we will 
see, the analogy between the Boltzmann factor and the plasma and the eigenvalue PDF for 
a certain random matrix ensemble requires that this constraint be relaxed. Nonetheless, we 
will find that up to terms which vanish as a Gaussian, the support of the eigenvalue PDF 
is still the spherical annulus. It should be mentioned that this analogy assumes a particular 
transformation of the eigenvalues, which start out as points in the complex plane. The mapping 
from a point z = x + iy in the complex plane, to a point {6, (p) on the sphere, is carried out by 
the stereographic projection 

a 

z = 2Re''*' tan - . (2.22) 

We know from e.g. [TSl eqns. (15.126), (15.127)] that then 



2R\u'v — uv\ 
Consequently, with z := z/{2R), 



cos -\Z — Z\ COS — . 

2' '2' 



dS 



1 



'l + \z\'^/AR 



2^2 



dxdy. 



N 



Uivir^'iur'^'' 



1=1 



N 

Y\ \ukVj - UjVkfdSi ■ ■ ■ dSN = Y\.{t 



15 |2 
\Zl\ 



l<j<k<N 



X 



n 



(l_L\^,\2)l3qN/2+2+l3{N-l)/2 11 I^J 
^ ' I 'I ^ l<j<k<N 



Zk\ 



1=1 

ir 1 ■ ■ ■ dfi^. 



+ \zi 



5, 1 2 



I3QN/2 



(2.23) 



We remark that the spherical annulus bounded between the azimuthal angles 6q and ^^r-e 



Making use of (2.22), together with (2.4) it follows that 



maps, under the stereographic projection (2.22), to a planar annulus with radii tq and r, 



3 Analogy with a random matrix ensemble 

Let A and B he N x N random matrices, with entries independently chosen as standard 
complex Gaussians. It was shown by Krishnapur fM] that the eigenvalue PDF of A~^B is, up 



to normalisation, given by the RHS of (2.23) with /3 = 2, g = Q = 0. In this section a more 
general random matrix realization of (2.23) will be given, applying for /3 = 2 and arbitrary qN , 
QN e Z>o. 

To achieve this, two results from random matrix theory must be combined. In relation to 
the first, with an n x M , n > M , standard complex Gaussian matrix a, set A = a^a to form a 
so-called complex Wishart matrix (see e.g. |il51 Ch. 3]). Let X be an M x iV, A^ > M, standard 
complex Gaussian matrix, then set Y = A~^^'^X. We know from [26j that, up to normalization, 
the element joint probability density function of Y is given by 

1 



det(I + ytr)r^+Jv■ 

In relation to the second of the results, suppose W is an M x N random matrix with element 
PDF of the form glWW''). Also, let ?7 be an M x M unitary random matrix chosen with Haar 
measure. Then we know from |12] that with N > M and up to normalization the PDF of 
G = Ui^WW^f^ is given by 

{deiG^Gf-^'g{G^G). 

Let us choose W in the second result according to Y as specified in the first. This shows 
that the element PDF of G = UiY^Yy/"^ is proportional to 

The explicit value of the proportionality constant can readily be calculated. 



Proposition 3.1 Let (3.1) when multiplied by 1/M he correctly normalized. Then we have 

For this to he well defined we require N > M and n > M. 

Proof. With C = GG"^ and the eigenvalues of C written {\j}j=i,...,m we know that 

{dG)=c W {\k-\fd\...\M. (3.3) 

l<3<k<M 



Here c is independent of the eigenvalues and {dG) denotes the product of differentials of the 
independent real and imaginary parts. To determine c, suppose temporarily that G is a standard 
complex normal random matrix so that it has PDF 



^_M2g-TrGtG 



11 e 






(3.4) 



Converting now to the corresponding measures on both sides using (|3.3|) then integrating shows 

(3.5) 



1 



TV C 



/•OO /"OO 

Jo Jo i^„-^i.^»,f 



l<j<k<M 



Evaluation of the integral (see e.g. [15], Prop. 4.7.3]) now gives 



TT 



nS'r(i+j)r(2 + j) 



(3.6) 



With c determined, we can proceed to evaluate J\f using an analogous strategy. Thus after 



multiplying (3.1 ) by l/A/" so that it is normalized from the analogue of (3.4) by introducing the 
eigenvalues of G^G. We then use (3.3) to convert that equation into an equality of measures. 



Integrating both sides, then changing variables Xj = tj/{l — tj) (j = 1, . . . , M) on the RHS we 

obtain 

»i /.I 



1 



c 



dti . 



M 



.N-M I 



tj) 



n-M 



n (^^ 



-ti 



dtM J_ J_ tj 

'0 "'O j=l l<j<k<M 

The multi-dimensional integral herein is a special case of the Selberg integral (see e.g. |10], [151 



Ch. 4]). It's evaluation as a product of gamma functions together with (3.6) gives (3.2) D 



We seek the eigenvalue PDF implied by the element PDF (3.1), normalized according to 



Proposition 3.1 Of course G as defined above (3.1) is non-Hermitian, and the eigenvalues will 



lie in the complex plane (for reviews of aspects of the rich mathematical physics associated 
with this setting see [22], [39], [33], [151 Ch. 15.]). We will see that the eigenvalue PDF can be 



identified with the RHS of (2.23) in the case (3 = 2, N = M and qN, QN e Z>o arbitrary. 



Proposition 3.2 Let G be an M x M matrix with element PDF (3.1) and normalized by (3.2). 
The corresponding eigenvalue PDF is given by 



M 



c 



nir 



1^ |2(Af-Af) 



n 



i=i 



_|_ U.|2^ra+Af-M+l 

^ 3\ I l<j<k<M 



\Zk- z 



']\ ' 



where 



C = M\7r 



M 



M-1 

n 

j=o 



T{N - M + 1 + j)T{n -M + l+j] 
r(n + A^-M + l) 



(3.7) 



(3.8) 



Proof. We follow [29j (see also [TJl Prop. 15.6.1]). The first step is to introduce the complex 
Schur decomposition by writing G = URU^ where U is an MxM unitary matrix and R = A+T, 
with A = diag(2;i, . . . , zm) the diagonal matrix of eigenvalues and T strictly upper triangular. 



8 



To make the decomposition unique, we must order the eigenvalues (for example, according to 
their modulus) and choose U from the right coset of the unitary group U[M] := U{M)/Ud{M), 
where Ud{M) denotes the set of diagonal M x M unitary matrices. The corresponding volume 
form is given by {U^dU). For later reference we note that (see e.g. [I5l eq. (3.23)]) 

Af(M-l)/2 

11^=0 r(j + i) 

We know that the change of variables formula from G to U and R is (see e.g. [151 Prop. 
15.1.1]) 



M 



{dG) 



n \zk-z,\''\[dx^dy, n dTJ.dTi.iU^dU), 

^<j<k<M j=l l<j<k<M 

where Zj = Xj + iyj and dTjk = dTJf, + idTji^. To obtain the eigenvalue PDF we must multiply 



this by (3.1), together with its normalization, and itegrate over U and T. Thus the eigenvalue 



PDF of G is equal to 



1 



{u^du) n 



M 



\2{N-M) 



i=i 



n \^k-zj\ 

l<j<k<M 



X 



det(l + WRY+^ 



n '^^kdT, 



jk- 



l<j<k<M 

Let Vm_i be an (m — 1) X 1 complex vector, and set 



'"m,p 



{d^r. 



m— ly 



TT 



^_^T{p-m + l) 



Also, after writing R = R^ to indicate the size of i?, set 

1 



T{p) 



(3.10) 



(3.11) 



''m,p\Z,\^ • • • 5 Zfn 



, TT dT'dT' 



Then we know from P^ eq. (15.138)] that 



c 



\Z\i ■ ■ ■ 1 z^ 



m,p 



This allows the final integral to be evaluated as 

CM-l,n+N-l 

JlT 
1=0 



-'m— l,p— ll^'^l; • • • ; Zj^i—l)- 



M-l 



Ha 



\zm-i 



2\n+N-M+l ■ 



Substituting in (3.10) and simplifying using ( 3.11[ ), (3.9) and (3.2) gives (3.7). In the normal- 
ization (3.8), the ordering on the eigenvalues has been relaxed. D 



Comparing (3.7) with the RHS of (2.23) we see that they agree if in the latter we set (3 = 2, 
N = M and 

QN = N -M, qN = n-M. 

In Figure 1 we show numerically generated eigenvalues corresponding to the choice Q = q = 1, 
stereographically projected onto the sphere. This illustrates the eigenvalue density being, to 
leading order, uniform within the spherical annulus, and zero outside. 




Figure 1: Stereograpically projected eigenvalues of matrices with element PDF (3.1 ) and eigen- 
value PDF (3.7) in the case n = N = 20, M = 10, repeated 1,000 times. The marked circles 



are the theoretical boundaries of support for n = A^ = 2M and M — > oo. The sphere has been 
scaled to have radius 1. 



4 Free energy 



Let us return now to the plasma interpretation of (3.7). A primary quantity of interest is then 
the large N form of the dimensionless free energy. 



/3F^ = -logZ;v(/3), 



where ^Ar(/3) is the partition function 

1/1 NAf/3/2 



ZnW) 



f±V\-^K.f[R^ '" 



N\\2RJ 



1=1 



deiivir^'iuif'^'' 



(4.1) 



JJ \ukVj - UjVkf . (4.2) 

^<j<k<N 



We know from studies relating the two-dimensional Coulomb gas to the Gaussian free field [31 
that the large N expansion of log Zn should be of the form 



logZjv(/3) ~ ApN + B^N'/^ + ^\ogN + --- . 



(4.3) 



Here —A/^ is the dimensionless free energy per particle, —Bjs is the dimensionless surface tension, 
and X denotes the Euler characteristic of the surface (explicitly x = 1 fo^' ^ disk, x = 2 for a 
sphere, x = for an annulus). 



The fact that the leading term in (4.3) is proportional to A^ follows from the proof of the 



existence of the thermodynamic limit for jellium by Lieb and Narhofer [36j. This term is a bulk 
quantity, and so is independent of the geometry. We know from [1] in the case of a disk that 
for /3 = 2 

(4.4) 



10 



_pb 

271 



2' 



That this is indeed independent of the geometry has been illustrated by exact calculation in 
the case of the sphere [U] , for example. Again from exact calculations in the case of the disk at 
/3 = 2, the exact form of B2 is known. It is expected to be dependent only on the length of the 
boundary, and exact calculation in the case of semi-periodic boundary conditions [7] illustrates 
this. In the case of soft wall boundary conditions, when the mobile particles are not confined to 
the region initially assumed in the computation of the Boltzmann factor, it has been observed 
in exact calculations [20j that B2 = 0. As we are interested in the case of soft wall boundary 
conditions, we thus expect that 

B2 = 0. (4.5) 



Hence the formula (4.3) predicts that for the plasma confined to the soft wall spherical annulus 

N 



logZjv(2) 



2 'OS ^ + 0(1). 



(4.6) 



Starting with (4.2), standard integration methods (see e.g. [13 §15.3]) verify (4.6), and further- 
more allow us to explicitly compute the term 0(1). 



Proposition 4.1 With (3 = 2, the asymptotic expansion of (4-2) for large N reads 

1 



logZ,v(2) 



T^°s2^ + i2^°^rTg 



1 , q 
— log 

12 ^1 + g 



O 



A^ 



(4.7) 



Proof. Recalling (|2.2|), simple manipulation of (|4.2|) in the case (3 = 2 gives 
Zf,(2) 



(±\ 



N 



m\2Rj 



,^2K^Y[R' 



X 



n 



Ui 



1=1 

Vk_ 



2tt 



n I dOi ( cos 



Q^-s2N-l+2qN r _ Q^-sl+2QN 



2) 



I sin 



Un 



Vk_ 
Uk 



(4.8) 



Making use of the Vandermonde determinant formula 






det[xj ^]j,fc=i,...,iv, 



the readily verified orthogonality 



R' 



2-K 



-\ {-\ =27rR^6^^n g{9)( tan-] sinOdO 



0^2n 



'0 Jo ^"^ ^"^ Jo 

valid for general g, and the Euler beta integral written in the form 



7r/2 



sin^)--^Hcos^)^^^^rf^^^^" + ')^(^+'^ 



T{a + b + 2) 



the integral (4.8) can be factorized into a product of one dimensional integrals with gamma 



function evaluations to give 

Zn{2) = (27ri?)^e-2^^ JJ 



N-l 






(4.9) 



11 



A formula more immediately suited for asymptotic analysis can be obtained by introducing 
the Barnes G-function. This satisfies the functional equation G{z + 1) = r{z)G{z), and can be 
given meaning for all complex z. In particular, it is known that for general a, 



N-l 



nr(i + / + 



a) 



1=0 



G{N + a + 1) 
G{a + 1) 



(see e.g. [151 eq. (4.183)]) allowing (4.9) to be rewritten 

N^-2K. 1 G(iV(l + g) + l)G(iV(l + g) + l) 



Zn{2) = {2TTRy^e 



(r(iv(i + g + g) + i)) 



N 



G{NQ + l)G{Nq+l] 



(4.10) 



In (4.10), using Stirling's formula for the gamma function, the known asymptotic formula for 



the Barnes G-function 



logG'(a; + l) ~ ^ 



3 T 1 / 1 \ 

logx- -x'^ + -log27r- — loga; + C'(-l) + 0(-j 



(see e.g. [HI eq. (14) ]) and recalling the explicit form (3.8), the stated expansion (4.7) then 
follows. n 



deduce the normalization (3.8). 



We remark that integrating both sides of (2.23) and using (4.9) is an alternative way to 



5 Correlation functions 

Throughout this section, we will work directly with the variables in the complex plane as implied 



by the eigenvalue problem, and are thus considering the RHS of (2.23) in the case (3 = 2. 
The general structure of the latter, being of the form 



N 



^<j<k<N 



h{r) :-- 



^2QN 



1=1 



^l _^^2-)(Q+g+l)7V+l 



(5.1) 



tells us that the fc-point correlation function has the determinantal form 

P(fe)(ri, . . . , Ffc) = det[i^(r^, r^)];,,T,=i,...,fc, 
where the so-called correlation kernel K is given by 



^(rA^.r^,) 



vr 



{h{r,)h{r,)Y/^H{r,r,e'^'^ 



ri-i 



1 ^ 



dr 



(5.2) 



(5.3) 



(here {r,6) are the polar coordinates of r). This follows from a simple calculation using the 
method of orthogonal polynomials (see e.g. [IHl Prop. 15.3.1]). 
We seek a form of K suitable for asymptotic analysis. 
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Proposition 5.1 Let 



and 



where 



h{r) : 

J{a, b; z) :- 

B{a,b) := 



1 



(1 _^ j.2YQ+q+l)N+l 

1 r t"-^ 



B{a,b) Jo (l + t)''+" 



dt. 



- t-i ^^ T{a)T{b) 



(5.4) 
(5.5) 
(5.6) 



,0 (l + t)''+" r(a + 6)' 

Furthermore, use |5.5| ) to define 

H{z) = (Q + g + l)Ar(l + z)(«+''+i)^-i( J(QiV, (g + l)iV; 2) - J((Q + l)Ar, qN; z)). (5.7) 

/n terms of the quantities h and H we have 



K{y^,y. 



1 ,7 



M' '■!) 



-0-,h 



-{h{r^)h{r,)Y'^H{r,r,d'^'--'^' 



(5.8) 



up to a factor which does not contribute to (5.2). 



Proof. Comparing (5.8) to (5.3), we see that the task is to find a summation formula, by 



way of an integral representation, of the summation defining H[z) in (5.3). Straightforward 

(5.9) 



working establishes that the latter satisfies the first order differential equation 

H'{z) + a{z)H{z) = b{z) 

where 

1 fQN 



a{z) 
b{z) 



,^ , -(g + l)iV + l 

1 + z \ z 

QNTi{Q + q + l)N + l] 



z^-'T{{Q + q + l)N + l) 



z{l + z)T{QN + l)r((g + 1)N) 1 + z T{N{Q + l))T{qN) 
According to the method of integrating factors, choosing I{z) such that 

I'{z) = I{z)a{z) 



(5.10) 
(5.11) 

(5.12) 



allows (5.9) to be written 



dz 



{I{z)H{z)) = b{z)I{z). 



Consequently H{z) can be expressed in terms of I{z) and b{z) according to 

1 



H{z) 



Hz) 



b{t)I{t)dt + C 



(5.13) 



for some C independent of z. 



Solving (5.12) gives 



yQN 



Hz) 



(l + ^)(Q+<?+i)Af-r 
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We substitute this and (5.11) into (5.13), then take the hmit z -^ to deduce that C = and 
thus conclude 



Hiz) 



(1 _^ ^\(Q+g+l)Af-l 



yQN 



QN 



^QN 



-1 



,dt 



B{QN + l,{q + l)N) Jo (l + t)(Q+«+i)^ 

qN f' t{Q+i)N~i 



B{{Q + l)N,qN+l) Jo (l + t)(Q+9+i)^ 



,dt 



(5.14) 



Use now of the recurrences 



B{x + l,y) 



X + y 



B{x,y), 



B{x,y + l) 



X + y 



B{x,y) 



in (5.14) gives the form (5.7), but with an extra fac tor o f l/z*^^. This latter factor is essentially 

ed by 

(5.15) 



cancelled by the factor of r^^ in {h{r)Y/'^ (recall (5.1)) in the sense that with h specified by 

H[z) = z^^H{z) 



(5.4) and 



we have that 



det 



Ch{r^)h{r^)Y/^H{r^r^e'^'^-'-'^ 



fj,,-y=l,...,k 



Since the above working shows that the formulas (5.15) and (5.7) for H(z) are consistent, we 
have established (5.8). D 



5.1 Global scaling 



In the variables of the RHS of (2.23), we know from (2.24) that the support of the underlying 



background charge density is between radii fg and fg, which are independent of N. Further- 
more the uniform background on the sphere maps, under the stereographic projection, to the 
background in the plane 

N{l + Q + q) 
-Pb[r) = _,, , nxo ■ (5.16) 



7r(l + r 



On the sphere, according to Proposition 2.1 the background density specified as the uniform 
value —ph within the spherical annulus and zero density outside is the solution of the integral 
equation 

— / p{{0, 4>)) log \uv — uv'\ dS = Cn — A^Qlogsin — — Nq cos —. 
Js 2 2 

As such ph provides the minimum of the energy functional 



E[P] 



a Qs 

p{{e, 0)) ( NQ log sin - + Nq log cos - ) dS 



dSip{{9i,(t)i)) / dS2p{{02,(t>2))^Og\u2Vi-UiV2\ 
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On the other hand, we know that to leading order the density of the mobile particles in the 
plasma can be characterised by minimizing this same energy functional (see e.g. |1] and ref- 
erences therein). Thus to leading order it must be that that the particle density is equal to 
Pb- When projected to the plane, this means that to leading order the particle density will 



be confined between radii tq and r^, and will have profile given by (5.16) (without the minus 



signs). We will see that this prediction is confirmed by explicit calculation, and we will show 
too that the correction terms are exponentially small in N . 



According to Proposition 5.1 



P(i)( 



{Q + q + l)N 



J{QN, {q + l)N- r^) - J{{Q + l)N, qN; r^ 



7r(l + r2)2 
Our task is to compute the large N asymptotic form of this expression. 



(5.17) 



Proposition 5.2 For asymptotically large values of N the density (5.11) vanishes outside the 
annulus r G [?"Q,Tg] up to exponentially small terms in N, while inside this annulus, again up 
to exponentially small terms P(i)(r) = Pbij) as specified by (5.1). 

Proof. According to (5.17) we require the large N form of J{aN, PN; x) for a; > fixed and 



a,/3 > 0. From the definition (5.5) we see that the A^-dependent portion of the integrand in 



the definition of J{aN, 13N; x) can be written 

N{alogt-{a+l3)log{l+t)) 



(5.18) 



This has a single maximum at t = to •= <^/P, and correspondingly J{aN, PN; x) is exponentially 
small when to is not part of the range of integration. Consequently, up to exponentially small 

terms in A^ 

1, X > a/ (3 

0, X < a/ (3. 

D 



J{aN, I3N; x) 



(5.19) 



The stated result now follows by using this result in (5.17). 



We now turn our attention to the large A^ behaviour of the truncated two-point correlation 
function, 

Pf2)(ri,r2) :=P(2)(ri,r2) - P(i)(ri)p(i)(r2). (5.20) 



According to (5.2) and (5.3) (for later purposes this is more useful than (5.8)) this has the 
explicit form 



Pf2)(ri,r2) 



TT^ 



;h{ri)h{r2)\H{rir2e 



ii9i-e2)\\2 



(5.21) 



For ri 7^ r2 and fixed as A^ — )■ oo, on the scale of the spacing between eigenvalues the eigenvalues 
at ri and r2 are effectively an infinite distance apart. They will thus be uncorrelated and so 
we expect p^Jri, r2) — ?■ as A^ — ?■ oo. On the other hand we now have that with ri = r2 the 
truncated two-particle correlation is equal to — (p(i)(ri))2 which we know is proportional to A^^ 
for ri inside the annulus. 

To quantify this behaviour, consideration of fiuctuation formulas for linear statistics (see 
Section 5]2] below) suggests that the appropriate quantity to analyze is 



/M:^/*,a(r0/dr,aM(p5,(r..r,)+.(r,-r,W,(r, 



for all a(r) sufficiently smooth. 
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Proposition 5.3 Let a(r) = a{r) = alr"^) so that a(r) is rotationally invariant, and with 



s, let a{s) be twice dijferentiable with respect to s We have 



lim I [a] 



(Q+l)/9 



{a'{s)Yds. 



(5.22) 



Q/a+q) 



Proof. We see from (5.3) that 

'■27r r'in 



dOi / d92\H{rir2e 
Jo 



i{9i-e2)\\2 ^ ^2 



N 



[rir2 



1 20-1) 



1 iir hirV'-^dry 



Consequently, with a(r) = a{r'^ 

I[a] = 



>^ //g°°Q;(r^)/i(r)r^-^ Mr 

i=i 



J h{r)r'^^~^dr 



v^ /Q°°(Q;(r^))^/i(r)r^-' Mr 



(5.23) 



Writing j = A^t, t := {j — 1)/N, and thus < t < 1, the large A^ form of the integrals can be 
determined as in the proof of Proposition 5.1 In particular, after changing variables s = r^, 
the maximum of the iV-dependent factor of the integrands (i.e. h{r)r^^~'^) is seen to occur at 



s = So(t) 



Q + t 
q + l-t' 



Thus we expand 



Mv^)4'"'^e-^(-^")'/2'^' 



a 



{Q + t){q + Q + l) 
(g+l-t)3 



\2„/// 



a[s) ~ a{so) + [s - so)a [so) + -{s - so) a [sq). 



Substituting in (5.23) allows us to conclude that for large N 

N 



^H ~ 2^(" (so)) ^^ 
7=1 Jo 



-7V(s-so)2/2aMs 



(g + Q + 1) ^, ,,2 JQ + t) 



But this last expression is just the Riemann sum approximation to an integral. After changing 
variables, (5.22) results. D 



We remark that for a{s) as in (5.3) 

"(Q+i)/g 

{a'{s)fds 

Q/(l+9) 



1 

2 
1 



''Q 



~~r a{r)dr 
r drJ 



['■Q'"''?] 






(5.24) 



16 



where Dy^Q^rq] denotes the annulus with inner radius tq and outer radius r^. This is consistent 
with the expected large A^ form ^30] 



Pf2)(ri, r2) + -^(ri - r2)p(i)(ri) = V?i^(ri - 1-2) 

for ri and V2 away from the boundary of the annulus. 

Coulomb gas theory predicts very different behaviour for ri, V2 within the boundary layer 
of the support pO]. Consider for definiteness the inner edge. The theory of [30] predicts 

where g{s\^ S2) has the property that 

/OO /"OO 

dsi I ds2g{si,S2) = 1. 

This result has previously been exhibited for the one-component plasma at /5 = 2 in the case 
of disk geometry [8], as has the analogue of (5.25) for the same system but now in an ellipse 
geometry [18] (the latter is equivalent to the partially symmetric Ginibre ensemble of complex 
random matrices [21]). It can readily be checked in the present setting of a projected spherical 
annulus. 



Proposition 5.4 The limit formula (5.25) holds true with 



9isuS2) = -e~''-'''K (5.26) 



n 



Proof. Our main tool is an asymptotic formula for J{aN, PN; x) valid for x bounded away 
from the real axis. In this case, along a ray from the origin to x, the corresponding integrand 
oscillates rapidly and the main contribution to the integral comes from the neighbourhood of 
the end point at x. To determine the latter we follow a strategy used on the incomplete gamma 
function in [8], involving a particular integration by parts. 

The integration by parts in turn is initiated by writing the integrand in terms of a derivative 
of its own functional form, 

r^-i 1 + 1 d / r^ 



(1 + t)("+/3)JV N{a - pt) dt\{l + t)(«+/3)^ 



Ul+tKa+/3)ivj- 



With X bounded away from the real axis, substitution of this in the definition (5.5) and inte- 
gration by parts shows 

1 1 x"^ / / 1 \ \ 

^'°^-^^'^'~ a(«iv.giV)iV(.-W(i + :^)"--^'" ^^''(^))- <"" 

Furthermore, with 7 := a//3, use of Stirling's formula shows 



Bi.N,0N)^J-{^) ( ) . (5.28) 



+ 7/ VI+7. 
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Substituting (5.28) in (5.27), then substituting the result in (5.14) and recaUing (5.21) shows 

N{Q + q + lf 



P{2)i.ZuZ2) 



\ZlZ2 



2QN 



27r3Q ((1 + |zi|2)(l + \z2\^))iQ+i+m+i 



l + rf.\ 2QN 



[l + r 



2\2iq+l)N-l 



1 - 2:1^2/^^ 



(5.29) 



Next, we must substitute for zi and Z2 as required by the LHS of (5.25). Appropriate large- iV 



expansion of the resulting terms on the RHS of (5.29) gives the RHS of (5.25) with g{si, S2) 
therein given by (5.26). D 



5.2 Local scaling 

A feature of the global scaling of the previous section is that the area of the annulus remains 
fixed as the number of eigenvalues tends to infinity. In contrast, we know from Section |4] that 
the thermodynamic limit is such that the volume of the annulus tends to infinity while the 
density of the eigenvalues stays fixed. At the level of the correlation functions, due to the scale 
invariance of the logarithmic potential, the thermodynamic limit is equivalent to a local scaling 
in which the position variables are measured on the scale of the (linear) inter-particle spacing. 
This can be achieved by rewriting each polar coordinate (r, 6) in terms of a cartesian coordinate 
{x, y) according to 

^ (5.30) 



X + 



X 



v^piPQ' 
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XV MX) 



kernel (5.3) under this scaling. 



for X e (rg, Tg), where Ph{,X) is given by (5.16). We seek the asymptotic form of the correlation 



Proposition 5.5 Let the polar coordinates of v^, r^ he replaced by the scaled cartesian coordi- 
nates (5.3(^. The correlation kernel (5.3) 



Ph{X) 



K{y^,,y^ 



,iNX{y^~y-,)/{^/p^{l+X'^)) 



TT 



TT 



X exp( - -(x^-a;^) --(?//.- 1/7) + «7r(a;^ + x^)(y^ - y^) + O 



1 



and thus, up to terms Oil/N^^/"^) 



Pt{X) 



Pik) 



>!,..., rfc) ~det[e 5(^m ^7)' 



^liyti-yiY+ii^iVti.-xti.y-t)' 



/i,7=l,...,fc- 



Ari/2. 
(5.31) 



(5.32) 



Proof. W e recall that Kivi^^, r^) is given in terms of h and H according to (5.8). But it follows 



from (5.19) that for r^, r^ in the annulus and within 0(1/ vA^) of the real axis 



(5.33) 



up to terms 0(1 vA^). Recalling now the definition (5.4) of h we thus have 



1 



Pb{X) 



K{r^,Y^ 



1 + r.rT,e*(''''"^^) \ (Q+q+i)N-i 



;i+r2)V2(l+ ^2)1/2 
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up to terms 0{l/yN). Introducing (5.30), the form (5.31) now follows upon elementary com- 



putation. And this used in (5.2), after observing that the first exponential factor on th e RHS 
does not contribute to the determinant, nor does the factors e*'^'^'''^''"^''^^''^ implies (5.32). D 



We would expect that the correlations in this bulk scaling limit would be independent of 
the geometry, and thus be the same as for the disk geometry for example. With z = x + iy the 
latter are given by [13 Prop. 15.3.2] 



detie-'^^l^^l +1^-1 )/2e^''^-]^,^=i,...,fc 



which is indeed identical to the RHS of (5.32). 



In (5.30) we required that X be strictly inside the annulus. With this assumption we were 



able to make use of (5.19). A physically different regime is to scale coordinates to have 0(1) 



spacing in the neighbourhood of a boundary of the annulus (for definiteness t his w ill be taken 
to be the inner boundary). With |zp = tq + 0{l/yN) the function J in (5.19) exhibits a 



crossover function form linking the two limiting values exhibited in (5.19) 



Proposition 5.6 Let J he specified by (5.5). We have 



lim J aN, (3N: - + , , 



1 1 c( 

- + -erf 

2 2 V 



^2) 



(5.34) 



where erf(a;) := -j^ f^ e * ds denotes the error function and 



Co 



a{a + (3) 



(5.35) 



Proof. With Co as in (5.35), we see that expanding the exponent on the RHS of (5.18) about 



its maximum at t = to to second order gives 



t 



aN 



laN 



-NCo{t-tof/2 



(l+t)(a+/3)JV (l+to)(«+m 

Recalling the definition (5.5[) of J, it follows that 



<'"''^''^'r^) 



X 



-'"/^dt 



which implies (5.34). 



Writing 



rQ + 



X 



+ 1- 



Y 



D 



(5.36) 



so that the (complex) coordinate is scaled and centred about the inner boundary, it follows 



from (5.34) that 



lim J{QN, (q + l)N; z^z. 



fi'-jj 



1 1 /x^^+x^+|(Vz^ 

2 2V ^ 



(5.37) 
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Proposition 5.7 Let ri, . . . ,rfc be centred and scaled about the inner boundary of the annulus 

e 

P(fc)(ri, ...,Vk)= det[H{{X^,Y^), iX^,Y^))U^=i,...,k (5.38) 



as implied by (5.36). We have 

1 \k 



lim 

TV^oo \pb[rQ] 



where 



HiiX,, Y,), (X„ Y,)) = e-|(^.-^.)^-|(^.->'.)^+*(^.^.-^M>'.) f 1 + ^erf f ^^ + ^^ + '^^^ ^^' 



2 2 



72 



(5.39) 



Proof. The only difference between this and the proof of Proposition 5.5 is that the function 



Hence the only difference between (5.38) and (5.32) is this extra factor 



J which implicitly multiplies the RHS of (5.33) is no longer unity but rather is given by (5.37) 



D 



The scaled edge correlation function (5.38) is precisely the same as found for the scaled 



edge correlation in the Ginibre ensemble of complex Gaussian matrices [ITj, which in turn 
is equivalent to the one-component plasma in an annulus at /3 = 2 with soft wall boundary 
conditions. 



5.3 Fluctuation formulas for linear statistics 

Knowledge of the one and two point correlation functions in the global scaling regime gives 
information on the mean and variance of a linear statistic. The latter is specified as the 
random variable A = Yli=i^i^j) where {zj} are the eigenvalues of the random matrix. Thus 
we have 

(A) = / p(^i){z)a{z) dxdy (5.40) 



and 



VarA 



dxidyi a{zi) / dx2dy2 a{z2) ( P(2)(^i, ^2) + <5(^i - Z2)p[i){z) 



(5.41) 



(for the latter equation, see e.g. [15, Prop. 14.3.2]), where pT^Jzi,Z2) is specified by (5.20). 
As previously remarked, the global scaling regime corresponds from a statistical mechanics 
viewpoint to an infinite density limit, since the N eigenvalues are confined to an annulus of 
fixed radius as A^ increases to infinity. It is precisely this limit (see e.g. [iSl Ch. 14]) that gives 
rise to universal behaviour by way of a Gaussian fluctuation formula, with a variance which is 
0(1). 

We will consider first the limiting form of the mean and variance. For the mean, it follows 



from (5.16) and (2.24) that 



lim —(A) 



[l + Q + q) 



TC 



a{z) 



^'-Q.-d 



(1 + kl 



2\2 



dxdy, 



(5.42) 
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where Z^j^p^fq] is defined as in (5.24). In the case that a{z) = a{\z\) so that the hnear statistic 
is rotationally invariant, we have from (5.22) with a(r^) = a{r) that the hmit variance is given 
by 

hmVarA=/ {a{yG)fds (5.43) 

^^°° JiQ+l)/q 

In this latter case the full distribution of A can easily be obtained via an explicit calculation. 
This is analogous to the situation for the complex Ginibre ensemble 



Proposition 5.8 Let (■) denote the average with respect to the PDF corresponding to the RHS 



of (2.23) in the case j3 = 2. Let a{z) = a{\z\) have a continuous derivative with respect to 
r = \z\ for r in D[f:q,fQ]- For large N we have that 



Qi!'T,iLia{zi) 



exp{zkN{q + Q + l) / 77^^^ 

{Q+i)/q U + ^) 

— lim Var A ) . 

2 Af-5-oo / 



p fQ/iq+l) 



(Q+i)/g 



a'(v^))2rfs 



exTpiikN lim — (A) 



(5.44) 



Consequently, as N ^ oo, A — (A) is distributed as a standard Gaussian with variance 
limiv_^oo VarA. 



Proof. Using an analogous integration procedure to that used in deriving (4.9), now using 



polar cordinates as in the workings of Sections |5.1| and 5.2, we readily obtain 



r"00 

N I e" 



Jf'T.t^io-i^i) 






VI 



1 



;; 



ds 



1 + s 



NQ / 1 NAr(g+l)+l 



I 



1 



s^-^ ds 



(5.45) 



We set j = Nt, where t := j/N and thus < t < 1. Now writing the A^ dependent terms 
in the integrand in the exponential form 



exp ( A^ 
we see that the maximum occurs for 



s 1 

Q log ^--- + (g + 1) log -—- + 1 log s 

1 + s 1 + s 



So(t) : = 



Q + t 
q + l-t' 



(5.46) 



(5.47) 



Expanding (5.46) to second order about So(t), and expanding the factor e'^^"'^^^ in the integrand 



of the numerator of (5.45) to first order about s = So{t) then completing the square, we see 
that 



Recalling now the definition of t above (5.46 ) we see that the sums are to leading order Riemann 



integrals. After a change of variables, the first line of (5.45) results. The second line follows by 
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using (|K42| and §^2^. D 



For linear statistics not rotationally invariant, there will be a contribution to the variance 



due to the universal form of the surface correlations (5.25) [HI [H]. In the case that a{r) is 
sufficiently smooth, a proof of its explicit form, together with a proof of the corresponding 
Gaussian fluctuation formula, follows from a more general theorem of Ameur, Hedenmalm and 
Makarov ^. This latter theorem also includes the setting of non-rotationally invariant linear 
statistics for the complex Ginibre ensemble, flrst established by Rider and Virag ^7\. In the 
case of linear statistics dependent only on the angle (and thus not smooth at the origin), the 
variance is typically no longer of order one, but nonetheless Gaussian fluctuation formulas can 
still be established [lOj . 
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